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Abstract. We derive an expression for the macroscopic force density that a narrow- 
band electromagnetic field imposes on a dissipative isotropic medium. The result is 
obtained by averaging the microscopic form for Lorentz force density. The derived 
expression allows us to calculate realistic electromagnetic forces in a wide range of 
materials that are described by complex-valued electric permittivity and magnetic 
permeability. The three-dimensional energy-momentum tensor in our expression 
reduces for lossless media to the so-called Helmholtz tensor that has not been 
contradicted in any experiment so far. The momentum density of the field does not 
coincide with any well-known expression, but for non-magnetic materials it matches 
the Abraham expression. 
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After the first well-known theoretical model for static electromagnetic forces in 
ponderable media, introduced by H. von Helmholtz [IJ, several other descriptions have 
been proposed to account for the action of both static and oscillating electromagnetic 
fields on a medium. The two most famous of them were proposed by H. Minkowski [2] 
and M. Abraham p]. Their models give different predictions for the value of the force 
density in certain particular cases, which has created a century-long scientific debate on 
the topic |1]-[15]. It is well known that both these models ignore the existing electro- 
and magnetostrictive forces [TJ HE] . A model by Einstein and Laub [I7j was intended to 
include these forces, but it turned out to be in disagreement with experiments (see, 
e.g., Ref. [7] and references therein). A reliable way to also include electro- and 
magnetostrictive forces can be based on the so-called Helmholtz energy-momentum 
tensor. Previously we have shown that the Helmholtz tensor can be derived in a 
relatively simple way starting with the microscopic field-matter interaction picture [H] . 
However, this tensor and the other models discussed so far are not applicable to 
dissipative media and, therefore, they do not allow calculating electromagnetic forces in 
many realistic materials. 

The aim of this work is to derive a general expression for the macroscopic force 
density imposed on a dissipative and electrically conducting inhomogeneous isotropic 
medium by a narrow-band electromagnetic field. We do this by spatially averaging 
the microscopic Lorentz force density. The result is written in terms of the three- 
dimensional electromagnetic energy-momentum tensor and the momentum density of 
the field, which explicitly depend on the complex relative permittivity and permeability 
of the material. In the limit of lossless medium, the obtained energy-momentum tensor 
becomes the Helmholtz tensor. We note that the force density expressed in terms of 
complex- valued quantities has been introduced previously (see, e.g., [19j and [20]), but 
not in connection with the Helmholtz tensor. The expression obtained by us for the 
field momentum density does not coincide with any of the well-known expressions, but 
for lossless dielectrics it converges to an expression obtained by us previously [18], and 
if the material is not magnetic, it matches the Abraham expression. 

The derivation of the force density in dissipative media is similar to that introduced 
by us in Ref. [18] for lossless materials. Since essential new details appear along the 
derivation, we present it here in its entirety. Within classical electrodynamics, the 
interaction of an electromagnetic field with matter is unambiguously described by the 
microscopic Maxwell's equations. 



eoV ■ e = e , 
V-b = 0, 
- V X e = b, 

— V X b = eoe + j , 



(1) 
(2) 
(3) 

(4) 



and the microscopic Lorentz force density 



mic 



^e+j X b. 



(5) 
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Here e and b are the microscopic electric and magnetic fields in the medium and e and b 
are their time derivatives. The electric charge and current densities ^ and j, respectively, 
can be written in terms of point electric charges qi as 

e = 5^g,5(r-r,), (6) 

i 

j = '^(liriS{r -Ti) , (7) 

i 

where 6{r — r^) is the Dirac delta function centered at the coordinate rj of charge qi. 

The bound electric charges in a medium can be combined into localized groups that 
belong to individual atoms (or molecules). For each such group one can expand the 
charge and current densities into Taylor series around the group's center r/ and then 
truncate the series to include contributions from electric and magnetic dipole moments 
d|^'' and m[^'' only [21]. Within this approximation the atoms are treated as point 
dipoles and the bound charge and current densities in the medium are 

^b = -Ed!^^-^'^(r-rO, (8) 

I 

Jb = E - + E ^ X - rz) , (9) 

I I 

The moments dj^'' and m|^^ contain the contributions of all bound electric charges of 
atom /. The second term in the expression for j^^ originates from the electric current 
loops due to rotational motion of the charges in the atoms [21] . 

If the medium contains conduction charges, an oscillating external field will put 
these charges into oscillating motion. For an electrically neutral medium, the oscillation 
of the conduction electrons is accompanied by an out-of-phase oscillation of the residual 
positive charges of the atoms, which can as well be treated as conduction charges. 

The oscillation of such positive and negative conduction charges produces their own 

(c) (c) 

dipole moments and . If the oscillating electron does not move far from 
a certain positively charged atom I (which is the case for high-frequency fields), the 
conduction charge and current densities can be written in the same form as ^-^ and j]^, 

i.e., ^c = -Ei ■ V5(r - n) and jc = Ei ^^^ir - r,) + E/ V x mf^6{r - r,). 

In order to take into account possible dissipation of electromagnetic energy in the 
medium, the moments d^^^ ni[^'', d['"'', and m[*^^ are considered to be complex vectors. 
For example, when calculating an expression for the frequency-dependent dielectric 
constant e, the contribution of each bound charge k to an atomic dipole moment turns 
out to be proportional to e~*'^*/('^fc — — ioj'jk), while the dipole moment originating 
from each conduction charge is proportional to e~*^*/(— tu^ — iowfj) |22]. Here and 
7fc are the resonance frequency and damping constant of the electron, respectively. For 
conduction electrons, Uj are equal to zero but jj are not because of ohmic losses. From 
this point of view the division of electric charges into bound and conduction charges is 
conditional. The total charge and current densities, ^ = + C,c and j = jjr, + jc are 
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Wc = o^^if miJ = T^Reire + j* x b}, (12) 



now given by the complex quantities 

e = -5^d,.V5(r-rO, (10) 
I 

j = J2 - rO + J] V X m^5(r - r^, (11) 

I I 

where d, = d^^^ + df^ and m, = mf^) + mf ^ are, respectively, the total electric and 
magnetic dipole moments per atom. 

When in addition to ^ and j the fields e and b are replaced with complex oscillating 
functions, the form of the microscopic Maxwell's equations and Eqs. (fTOj) and (ITTIl 
are preserved, but the Lorentz force density must be rewritten. In terms of complex 
quantities, one can write the microscopic Lorentz force density averaged over one 
oscillation period as 

-Re{/ j2iic} = 2' 

where Re denotes the real part, and the asterisk stands for complex conjugation. Here 
the functions e and b are assumed to oscillate harmonically with slowly varying complex 
amplitudes, in which case Au <^ u, with Aoj and oj denoting the field's bandwidth and 
carrier frequency, respectively. In this case, the changes in the field amplitudes during a 
single oscillation period of the field are negligibly small. In Eq. ( lT2l) the fast oscillation 
at the carrier frequency is averaged out and fj^j,^ remains a slowly varying function of 
time. 

Using Eqs. (fTOj) and (fTTI) . the function /j^i^ is written as 

/mic = E (- ■ "^^(^ - ^')^(^) + - rz) x b(r) (13) 

+ (V X m^6{r — r;)) x b(r) 

By integrating this function over a representative elementary volume 6V and dividing the 
result by this volume, one can find the macroscopic force density function / = (/jxiic); 
where the angle brackets denote spatial averaging. The function / is given by |18] 

/ = ^ E (v(d;.e(r,)) + V(m;.b(rO) + ^(drxb(rO)y (14) 
I in sv 

where the summation is performed over atoms belonging to 6V. 

Equation (fT2ll implies that the macroscopic force density averaged over one 
oscillation period of the field is 

f=(fmic) = ^M/}- (15) 

Note that the order of temporal and spatial averaging does not affect the result. For 
charges within 6V and for an arbitrary coordinate r in 6V, one can write e(r) = 
eext(r) +eint(r) and b(r) = bext(r) + bint(r), where e^^^{r) and h^^^{r) are the fields 
produced by the charges in SV and the fields egj^^(r) and bgxt(r) are created by all 
sources that are external to 6V. The external fields are independent of the charges and 



Electromagnetic force density in dissipative isotropic media 



5 



their coordinates in 6V, while the "internal" fields are strongly inhomogeneous around 
the internal charges. Substituting these expansions into Eq. (|T^ and using Eq. f fTSj) . 
one obtains 

f=^Re{^ E (v(d;-eext(rO) + V(mr-bext(rO) (16) 
I in SV 

+ X bext(rO))} +fint' 

where the second term, i^^f^, has the same form as the first one, but with the external 
fields replaced with the internal ones [18]. The term fj^^^ is equal to 1/6V times the 
total force imposed on all electric charges in 6V by the spiking fields ej^^ (r) and bj^^^ (r) 
produced by the charges themselves. According to Newton's third law, this force is 
equal to zero and thus fj^^^ = 0. Removing fjj^|- from Eq. f lTB]) and assuming that in the 
small volume 6V the total electric and magnetic dipole moments of each atom are equal 
to constant vectors d and m, one obtains an expression for f = Re{/}/2, where / is 
given by 

/ = ^E(^^^ E ^ext..(rO+<V E W(rO) (17) 

k=x,y,z I'mSV I in 5V 

+ ;|(d'x^ E bextW). 

I in SV 

Here the quantities with subindex k are the Cartesian components of the corresponding 
vector quantities. 

The atoms can to the first order be considered to be distributed uniformly within 
the small 5V . Since in 5V the fields eg^t (r) and bg^t (i") ^'^^ smooth charge-independent 
functions, one can replace the averaging of the fields over the coordinates of the atoms 
in Eq. ( fT7|) with volume averaging, which results in 

f=^^H 5Z (^;^*VEext,. + ^:V5gxt,fc)+^(P*xBext)}, (18) 

k=x,y,z 

where Eg^t = (Gext('^)) ^^"^ -^ext = (bext('^))- '^^^ electric polarization P and 
magnetization M are given by 

V = Nsvdl6V = V^ + Vc. (19) 
M = Nsvm/5V = + Mc , (20) 

with N^v denoting the number of atoms in 5V . In the above equations, P]^ and 
originate from bound charges and Pc and Mc from conduction charges. We note that 
Eq. f|T8l) holds also for anisotropic and nonlinear materials. In what follows, however, it 
is assumed that the medium is isotropic and linear, so that the polarization components 
are described by P]q = eo(ei3 — 1)E and Pc = 2Ecr/(eow), where is the complex 
dielectric constant due to bound charges and a the complex electric conductivity; i 
is the imaginary unity. Usually the polarization is written as P = eo(e — 1)E, where 
e = + ia/{eou) is the overall dielectric constant [22]. Similarly, the magnetization 
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can be written as M = (/i — 1)H, with /i being the overall complex relative permeability 
of the medium. The above equations are valid, if the size 5V^^^ is much larger than 
the distance le over which a conduction electron moves during one oscillation period of 
the field. This condition limits the theory to high-frequency fields. However, even for 
copper in an electric field with an amplitude of 1 kV/cm, the distance le is on the order 
of 1 i^m or less already at a frequency of 1 GHz. 

The averaged external electric field can be found from Eg^t = E — (ej^^^-) [18]. For 
an arbitrary charge distribution in 6V, the field (ej^^^-) is equal to —Dsv/i'i^oSV), where 
Dgy = P6V is the total dipole moment of the medium within SV (see also sect. 2.13 
in [23]). Therefore, the external field is 

Eext = E + P/(3eo) = E(e + 2)/3, (21) 

which has the same form as the traditional local field with the Lorentz correction. Note, 
however, that the polarization P also contains the contribution from the electric dipole 
moments of the conduction charges. Similarly, the external magnetic field is calculated 
as Bgxt = B — (bjj^^). The field (bj^^) created by microscopic electric current loops in 
6V is given by (bj^^) = 2/ioM/3 (see Eq. (9-22) in [23]), which yields 

Bgxt = B - 2/ioM/3 = /xoH(/i + 2)/3, (22) 

where the relation B = /io/iH has been used. Substituting the calculated Eg^t and 
Bgxt into Eq. (fTSl) and expressing P and M through E and H, one obtains the force 
density 

i=l^^{ ''^'*~ ^\ \E\'Ve+{e + 2) ^IVE,] (23) 



k=x,y,z 



k=x,y,z 



that depends on E and H and on the complex parameters e and /i. 

The macroscopic fields E and H satisfy the macroscopic Maxwell's equations 

V ■ D = , (24) 
V-B = 0, (25) 
- V X E = B , (26) 

V X H = D , (27) 

where the electric charge and current densities due to conduction charges are absorbed 
in D = eocE. With the help of these equations and using standard differential identities 
involving vectors and dyads, the obtained equation for f can be rewritten in a general 
form as 

^ HG 

f=-V-T-^, (28) 
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where the energy-momentum tensor T and the momentum density G are given by 



Here, E*E and H*H denote the outer products of the vectors, E and H are the complex 
amphtudes of the fields, i denotes the unit tensor, and c is the speed of light in vacuum. 
The form of Eq. fl28|) is not only physically insightful, but also very convenient in view 
of calculations, since for stationary fields the force on an object in a medium can be 
calculated simply by integrating n ■ T over the surface inclosing the object instead of 
integrating f over the volume of the object; n is the unit vector normal to the surface. 
It can be seen that if the medium is lossless, so that both e and /i are real, Eqs. ( 129|1 
and ( 130|) converge to Eqs. (35) and (36) in [18]. The tensor in Eq. (35) of Ref. [18] is 
the Helmholtz tensor (see, e.g., Ref. [?]), and the field momentum density of Ref. [18] 
still has to be experimentally verified. At high frequencies, n is equal to 1 for most 
materials, and the momentum density G becomes 



independently of e. Equation ( l3T]) is seen to match the Abraham's expression for G 
averaged over one oscillation period of the field within the slowly varying envelope 
approximation. Since Eqs. f l28|) -f p0|) are valid for dissipative media, they can be applied 
also to metamaterials that can have /i 7^ 1 even at optical frequencies. The quantity 
dG/dt is equal to zero for constant-amplitude harmonic fields, and in the case of a 
slowly varying amplitude the contribution of this quantity to f is small compared to the 
contribution of V ■ T. 

It is worth mentioning that if one deals with dielectric liquids at hydrodynamic 
equilibrium with the field, the negative of the gradient of the hydrostatic excess pressure 
can be found to compensate for the electromagnetic strictive force density at each point 
in the liquid. Inclusion of this hydrostatic force density often leads to a result that 
can as well be obtained by using Minkowksi's or Abraham's tensors [7]. However, in a 
non-equilibrium case, e.g., immediately after introducing the field in the medium, the 
pressure evolves and does not compensate for the electromagnetic strictive force. 

As a simple example of new phenomena that can be revealed by applying Eqs. (I28l) - 
( !30|) . let us consider a monochromatic electromagnetic plane wave (dG/dt = 0) that 
propagates along z-axis in a dissipative medium and has an attenuation length of zg,. 
At 2; = 0, the complex amplitudes of the field are Eq and Hq. The time- averaged force 
imposed by the field on the part of the medium that is confined between the planes 
z = and z = z' is 




(29) 



(30) 




(31) 




(32) 
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where the surface integrals are obtained by applying Gauss' integration law to the 
original volume integral; rii and 112 are the normals to the integration surfaces directed 
outward from the section of the medium in question. If E and H are everywhere 
perpendicular to z and 2;' ^ 2;^, Eq. (132|) becomes 



2 + 2Re{e} - |e| 
6 



\Eo\'^dxdy (33) 



6 

where Eq. f l29|) has been used; z is the unit vector along z. The integrals in 
Eq. f p3|) are proportional to the intensity Jg of the plane wave at z = 0, since 
l-^oP = 2/o^/io/eo|/i|/Re{^/e7i} and |ifoP = 2/oA/eo//io|e|/Re{^e/I}. In terms of Iq, 
the force F per unit cross-sectional area of the wave, pressure p, reads 

z/o 2 + 2Re{e}- |e|2 + 3|e| 

P = cR^ 6 ' 

where the medium is assumed to be non-magnetic so that fi = 1. 

It can be seen that for sufficiently large |e|, the pressure p becomes directed opposite 
to the propagation direction of the beam. Thus, if one would consider the matter-field 
interaction purely in terms of momentum exchange between photons in the wave and 
the medium, the obtained average momentum per photon would be negative. Obviously, 
photons can not only share their momenta with the medium, but also impose gradient 
forces on it. The pressure p in Eq. flM|) contains both a positive component due 
to the momentum exchange (the usual radiation-pressure component) and a negative 
component due to the intensity gradient originating from the attenuation of the wave. 
For large |e|, the second contribution can be stronger than the first one. Single-crystal 
silicon can be considered as an example. For a wavelength of, say, 532 nm (frequency- 
doubled Nd:YAG laser), we have e ^ 20.3 + il.O, Re{^/e} ^ 4.5, and |e| ^ 20.3 |2l]. 
Substituting these values into Eq. ( IMI) yields p = — 11.4z/o/c. This pressure points 
against z and is 11.4 times higher than a pressure that would result from full absorption 
of the same wave by an object in vacuum. Another example is salt water at 1-THz 
field frequency. For a mass fraction of 0.25 of NaCl in the solution, the solution is 
characterized by e ^ 4.9 + i3.1, Rej^e} ~ 2.3, and |e| ~ 5.8 [25]. The resulting pressure 
is p = — 0.3z/o/c that is also opposite to z. 

Suppose now that an electromagnetic beam instead of a plane wave is interacting 
with the previous medium of salt water. If the beam diameter is large compared to Za, 
then switching on the beam will lead to a motion of the liquid toward the beam source. In 
time, the hydrostatic pressure will be redistributed to compensate for the gradient force 
density due to the field. However, the positive radiation-pressure force component will 
remain uncompensated and, hence, the liquid will eventually fiow along the beam axis 
in the positive direction of z. In principle, the same motion of the hydrodynamically 
equilibrated liquid can be obtained also within Minkowski's and Abraham's theories 
written in the domain of complex functions, but what cannot be obtained using these 
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theories is the intensity dependent hydrostatic pressure. This is explained by the fact 
that Minkowski's and Abraham's force densities imphcitly contain the action of the 
medium on itself in the form of an Archimedes-like force density that at equilibrium 
compensates for the compressive action of the field. 

In conclusion, we have derived the equation for the force density imposed by a 
narrow-band electromagnetic field on a medium that is characterized by complex-valued 
electric permittivity and magnetic permeability. Equation fl28|) expresses the force 
density in the form that allows one to conveniently calculate the overall force imposed 
by the field on an arbitrary part of the medium by evaluating a surface instead of a 
volume integral, if the field is stationary. For narrow-band fields, Eq. ( |30l) can be used 
to evaluate the momentum density of the field averaged over its single oscillation pediod. 
This quantity depends on e and /i, but it takes on the Abraham's form if the material 
is not magnetic, which is the case for essentially all materials in high-frequency fields. 
However, Eqs. (!28|) -( l30|) can be applied also to high-frequency magnetic metamaterials, 
which can reveal unexpected phenomena associated with electromagnetic forces. 
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